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The Dirac cone underlies many unique electronic properties of graphen# and topological 
insulator^, and its band structure—two conical bands touching at a single point—has also 
been realized for photons in waveguide arrayiP, atoms in optical latticed, and through acci¬ 
dental degeneracy SEL Deformations of the Dirac cone often reveal intriguing properties; an 
example is the quantum Hall effect, where a constant magnetic field breaks the Dirac cone 
into isolated Landau levelP. A seemingly unrelated phenomenon is the exceptional point* 11 , 
also known as the parity-time symmetry breaking point 12 ls , where two resonances coincide 
in both their positions and widths. Exceptional points lead to counter-intuitive phenomena 
such as loss-induced transparency®, unidirectional transmission or reflection 17 23 , and lasers 
with reversed pump dependence®® or single-mode operation 27 28 . These two fields of re¬ 
search are in fact connected: here we discover the ability of a Dirac cone to evolve into a ring 
of exceptional points, which we call an “exceptional ring.” We experimentally demonstrate 
this concept in a photonic crystal slab. Angle-resolved reflection measurements of the pho¬ 
tonic crystal slab reveal that the peaks of reflectivity follow the conical band structure of a 
Dirac cone from accidental degeneracy, whereas the complex eigenvalues of the system are 
deformed into a two-dimensional flat band enclosed by an exceptional ring. This deforma¬ 
tion arises from the dissimilar radiation rates of dipole and quadrupole resonances, which 
play a role analogous to the loss and gain in parity-time symmetric systems. Our results in¬ 
dicate that the radiation that exists in any open system can fundamentally alter its physical 

properties in ways previously expected only in the presence of material loss and gain. 

Closed and lossless physical systems are described by Hermitian operators, which guaran¬ 
tee realness of the eigenvalues and a complete set of eigenfunctions that are orthogonal to each 
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other. On the other hand, systems with open boundaries 10 29 or with material loss and gain 12 14 EM 28 


are non-HermitiarP and have non-orthogonal eigenfunctions with complex eigenvalues where the 
imaginary part corresponds to decay or growth. The most drastic difference between Hermitian 
and non-Hermitian systems is that the latter exhibit exceptional points (EPs) where both the real 
and the imaginary parts of the eigenvalues coalesce. At an EP, two (or more) eigenfunctions col¬ 
lapse into one so the eigenspace no longer forms a complete basis, and this eigenfunction becomes 
orthogonal to itself under the unconjugated inner product 8 11 . To date, most studies of EP and 
its intriguing consequences concern parity-time symmetric systems that rely on material loss and 
gain 12 14 OH" 28 , but EP is a general property that requires only non-Hermiticity. Here, we show the 
existence of EPs in a photonic crystal slab with negligible absorption loss and no artificial gain. 
When a Dirac-cone system has dissimilar radiation rates, the band structure is altered abruptly 
to show branching features with a ring of EPs. We provide a complete picture from analytic 
model and numerical simulation to experimental observation; together, they illustrate the role of 
radiation-induced non-Hermiticity that bridges the study of EPs and the study of Dirac cones. 

We start by showing that non-Hermiticity from radiation can deform an accidental Dirac 
point into a ring of EPs. First, consider a 2D photonic crystal (PhC) 30 (inset of Fig. la), where 
a square lattice (periodicity a) of circular air holes (radius r) is introduced in a dielectric mate¬ 
rial. This is a Hermitian system, as there is no material gain or loss and no open boundary for 
radiation. By tuning a system parameter (for example, r), one can achieve accidental degeneracy 
between a quadrupole mode and two degenerate dipole modes at the T point (center of the Bril- 
louin zone), leading to a linear Dirac dispersion due to the anti-crossing between two bands with 
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the same symmetry 5 —. The accidental Dirac dispersion from the effective Hamiltonian model 
(see equation ([l} below with 70 = 0) is shown as solid lines in Fig. la , agreeing with numeri¬ 
cal simulation results (symbols in Fig. la). In the effective Hamiltonian we do not consider the 
dispersionless third band (gray line) due to symmetry arguments (section I in Supplementary In¬ 
formation), although this third band cannot be neglected in certain calculations, including Berry 
phase and effective medium property 32 33 . 

Next, we consider a similar, but open, system: a PhC slab (inset of Fig. lb) with finite thick¬ 
ness h. With the open boundary, modes within the radiation continuum become resonances because 
they radiate by coupling to extended plane waves in the surrounding medium. Non-Hermitian per¬ 
turbations need to be included in the Hamiltonian to account for the radiation loss. To the leading 
order, radiation of the dipole mode can be described by adding an imaginary part —i'fd t0 the 
Hamiltonian, while the quadrupole mode does not radiate due to its symmetry mismatch with the 
plane waves 34 . Specifically, at the T point the system has 62 rotational symmetry (invariant under 
180° rotation around the z axis), and the quadrupole mode does not couple to the radiating plane 
wave because the former is even [E(r) = Oc 2 E(r)] whereas the latter is odd [E(r) = —Oc 2 E(r)] 
under C 2 rotation 30 . The effective Hamiltonian is 


H cS = 


, \ 

UJ 0 Vgk 

Vgk UJ 0 - *7d 


( 1 ) 


with complex eigenvalues 


w± = LUo - ± V g ^/k 2 - k%, (2) 

where ujq is the frequency at accidental degeneracy, v g is the group velocity of the linear Dirac 
dispersion in the absence of radiation, k is the magnitude of the in-plane wavevector (k x , k y ), and 
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Ay = r )&/2v g . Here, one of the three bands is decoupled from the other two and is not included 
in equation (|T) (see section II of Supplementary Information). In equation ([2]), a ring defined by 
k = k c separates the k space into two regions: inside the ring (k < k c ), Re(u;±) are dispersionless 
and degenerate; outside the ring (k > Ay), Im(o,' ± ) are dispersionless and degenerate. In the vicinity 
of k c , Im(o;±) and Re(w±) exhibit square-root dispersion (also known as branching behavior) inside 
and outside the ring, respectively. Exactly on the ring (k = Ay), the two eigenvalues ui± are 
degenerate in both real and imaginary parts; meanwhile, the matrix H e g becomes defective with an 
incomplete eigenspace spanned by only one eigenvector (1, -i) T that is orthogonal to itself under 
the unconjugated inner product. This self-orthogonality is the definition of EPs; hence, here we 
have not just one EP, but a continuous ring of EPs. We call it an exceptional ring. 

Fig. lb,c show the complex eigenvalues of the PhC slab structure calculated numerically 
(symbols), which closely follow the analytic model of equation ([2]) shown as solid lines in the 
figure. When the radius r of the holes is tuned away from accidental degeneracy, the exceptional 
ring and the associated branching behavior disappear, as shown in Fig. SI. Several properties 
of the PhC slab contribute to the existence of this exceptional ring. Due to periodicity, one can 
probe the dispersion from two degrees of freedom, k x and k y , in just one structure. The open 
boundary provides radiation loss, and the C 2 rotational symmetry differentiates the radiation loss 
of the dipole mode and of the quadrupole mode. 

We can rigorously show that the exceptional ring exists in realistic PhC slabs, not just in 
the effective Hamiltonian model. Our proof is based on the unique topological property of EPs: 
when the system parameters evolve adiabatically along a loop encircling an EP, the two eigenvalues 
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switch their positions when the system returns to its initial parameters 10 11 29 35 , in contrast to the 
typical case where the two eigenvalues return to themselves. Using this property, we numerically 
show, in Fig. S2 and section III of Supplementary Information, that the complex eigenvalues always 
switch their positions along every direction in the k space, and therefore prove the existence of this 
exceptional ring. As opposed to the simplified effective Hamiltonian model, in a real PhC slab, the 
EP may exist at a slightly different magnitude of k and for a slightly different hole radius r along 
different directions in the momentum space, but this variation is small and negligible in practice 
(section IV of Supplementary Information). 

To demonstrate the existence of the exceptional ring in such a system, we fabricate large-area 
periodic patterns in a SisN 4 slab (n = 2.02, thickness 180 nm) on top of 6 /im of silica (n = 1.46) 
using interference photolithography' 4 . Scanning electron microscope (SEM) images of the sample 
are shown in Fig. 2a, featuring a square lattice (periodicity a = 336 nm) of air cylindrical holes 
with radius of 109 nm. We immerse the structure into an optical liquid and tune the refractive 
index of the liquid; accidental degeneracy in the Hermitian part is achieved when the liquid index 
is selected to be n = 1.48. We perform angle-resolved reflectivity measurements (setup shown in 
Fig. 2b) between 0 and 2 degrees along the T to X direction and the F to M direction, for both s 
and p polarizations. The measured reflectivity for the relevant polarization is plotted in the upper 
panel of Fig. 2c, showing good agreement with numerical simulation results (lower panel), with 
differences coming from scattering of disorder, inhomogeneous broadening, and the uncertainty in 
the measurements of system parameters. The complete experimental result for both polarizations 
is shown in Fig. S3; the third and dispersionless band shows up in the other polarization, decoupled 
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from the two bands of interest. 


The peaks of reflectivity (dark red color in Fig. 2c) follow the linear Dirac dispersion; this 
feature disappears for structures with different radii that do not reach accidental degeneracy (exper¬ 
imental results in Fig. S4). To understand the reflection peaks, we consider a generic two-by-two 
Hamiltonian H with no assumption made about its matrix elements. We separate H into a Hermi- 
tian part A and an anti-Hermitian part -zB (so that A and B are both Hermitian), and choose the 


basis in which A is diagonal: 

( 

uhu t = 
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(3) 


A /B 

As before, we use u>± to denote the complex eigenvalues of the Hamiltonian A — zB. The reflec¬ 
tivity in our system can be modeled using temporal coupled-mode theory (TCMT, with details in 
section V of Supplementary Information), where we show that the reflection peaks generally occur 
near the eigenvalues fi 1)2 of the Hermitian part A and are independent of the non-Hermitian part 
—zB (Fig. S5 with details in section VI in Supplementary Information). Therefore, the peak loca¬ 
tions in Fig. 2c (dark red) reveal information about only the Hermitian part of the Hamiltonian; the 
fact that they show linear Dirac dispersion indicates that we have successfully achieved accidental 
degeneracy in the eigenvalues of the Hermitian part, consistent with the simplified model in equa¬ 
tion O In Fig. S6, we plot the fi li2 extracted from the reflectivity data through a more rigorous 
data analysis using TCMT (described below); the linear dispersion is indeed observed. 

The eigenvalues of the Hamiltonian, u>±, behave very differently from the reflectivity peaks. 
Simulation results (white lines in the lower panel of Fig. 2c) show Re(u;±) are dispersionless at 
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small angles with a branch-point singularity around 0.31°—consistent with the feature predicted 
by the simplified Hamiltonian in equation [2] In Fig. 2d, we compare the reflectivity spectra from 
simulations (with peaks indicated in red arrows) with the corresponding complex eigenvalues at 
three representative angles (0.8° in blue, 0.31° in green, and 0.1° in magenta). At 0.31°, the two 
complex eigenvalues are degenerate, indicating an EP; however, the two reflection peaks do not 
coincide since they represent the eigenvalues of only the Hermitian part of the Hamiltonian, which 
does not have degeneracy here. The dip in reflectivity between the two peaks (marked as black 
arrows in Figs. 2 and 3) is the coupled-resonator-induced transparency (CRIT) that arises from the 
interference between radiation of the two resonance^ 37 , similar to electromagnetically induced 
transparency (HIT; 3 *. 

To extract the underlying Hamiltonian matrix and its eigenvalues from the measured re¬ 
flectivity spectrum, we use TCMT to model the direct and the resonant reflection processes; the 


expression for reflectivity is given in equation ( |S. 15[ ) with the full derivation given in section V 
of the Supplementary Information. Fitting the reflectivity curves with the TCMT expression gives 
us the matrix elements of the Hamiltonian (as shown in equation ([3])) that we use to calculate its 
eigenvalues; this procedure is the same as our approach in Ref. [39]except that here we handle mul¬ 
tiple resonances simultaneously, accounting for their non-orthogonality and radiative coupling 40 . 
Fig. 3a compares the fitted and the measured reflectivity curves at three representative angles (with 
more comparison in Fig. S6a); the excellent agreement shows the validity of the TCMT equations. 
Underneath the reflectivity curves, we show the complex eigenvalues. 

Repeating the fitting procedure for reflectivity spectrum measured at different angles, we 
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obtain the dispersion curves for all complex eigenvalues, which are plotted in Fig. 3b. Along 
both directions in k space (r —* X and T —> M), the two bands of interest (shown in blue and 
red) exhibit the EP behavior predicted in equation (|2|): for k < k c the real parts are degenerate 
and dispersionless; for k > k c the imaginary parts are degenerate and dispersionless; for k in the 
vincinity of k c branching features are observed in the real or imaginary part. In Fig. 3c, we plot the 
eigenvalues on the complex plane for both the r —>■ X and T —>• M directions. We can see that in 
both directions, the two eigenvalues approach each other and become very close at certain k point, 
which is a clear signature of the system being very near EP. 

We have shown that non-Hermiticity arising from radiation can significantly alter fundamen¬ 
tal properties of the system including the band structures and density of states; this effect becomes 
most prominent near EPs. The PhC slab described here provides a simple-to-realize platform for 
studying the influence of EPs on light-matter interaction, such as for single particle detection 41 
and modulation of quantum noise? 42 . The two-dimensional flat band also provides high density 
of states and therefore high Purcell factors. The strong dispersion of loss in the vicinity of the T 
point can improve the performance of large-area single-mode PhC lasers^ 43 . The deformation into 
exceptional ring can also occur for non-accidental Dirac point^S. Further studies can advance the 
understanding of the connection between the topological property of Dirac points 2 45 and that of 
EPs 35 in general non-Hermitian wave systems, and this method of our study goes beyond photonics 
to phonons, electrons, and atoms. 

METHODS SUMMARY 

Sample fabrication. The Si 3 N 4 layer was grown with the low-pressure chemical vapor deposition 
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method on a 6/rm-thick cladding of Si0 2 on the backbone of a silicon wafer (LioniX). Before 
exposure, the wafer was coated with a layer of polymer as anti-reflection coating, a thin layer of 
Si0 2 as an intermediate layer for etching, and a layer of negative photoresist for exposure. The 
square lattice pattern was created with Mach—Zehnder interference lithography using a 325-nm 
He/Cd laser. The angle between the two arms of the laser beam was chosen for a periodicity of 
336 nm. After exposures, the pattern in the photoresist was transferred to Si 3 N 4 by reactive-ion 
etching. 

Experimental details. The source was a supercontinuum laser from NKT Photonics (SuperK- 
Compact). A polarizer selected s- or p-polarized light. The sample was immersed in a colorless 
liquid with tunable refractive indices (Cargille Labs). The sample was mounted on two perpen¬ 
dicular motorized rotation stages (Newport): one to orient the PhC to the T-X or T-M direction, 
and the other to determine the incident angle 6. The reflectivity spectra were measured with a 
spectrometer with spectral resolution of 0.02 nm (HR4000; Ocean Optics). 


1. Castro Neto, A. H., Guinea, F., Peres, N. M. R., Novoselov, K. S. & Geim, A. K. The electronic 
properties of graphene. Rev. Mod. Phys. 81 , 109-162 (2009). 

2. Hasan, M. Z. & Kane, C. L. Colloquium: topological insulators. Rev. Mod. Phys. 82 , 3045 
( 2010 ). 

3. Rechtsman, M. C. et al. Strain-induced pseudomagnetic field and photonic landau levels in 
dielectric structures. Nature Photonics 7 , 153-158 (2013). 


10 




4. Tarruell, L., Greif, D., Uehlinger, T., Jotzu, G. & Esslinger, T. Creating, moving and merging 
dirac points with a fermi gas in a tunable honeycomb lattice. Nature 483 , 302-305 (2012). 

5. Huang, X., Lai, Y., Hang, Z. H., Zheng, H. & Chan, C. Dirac cones induced by accidental 
degeneracy in photonic crystals and zero-refractive-index materials. Nature materials 10 , 
582-586 (2011). 

6 . Moitra, P. et al. Realization of an all-dielectric zero-index optical metamaterial. Nature Pho¬ 
tonics 7 , 791-795 (2013). 

7. Zhang, Y., Tan, Y.-W., Stormer, H. L. & Kim, P. Experimental observation of the quantum 
Hall effect and Berry’s phase in graphene. Nature 438 , 201-204 (2005). 

8 . Moiseyev, N. Non-Hermitian quantum mechanics (Cambridge University Press, 2011). 

9. Kato, T. Perturbation theory for linear operators, vol. 132 (Springer Science & Business 
Media, 1976). 

10. Rotter, I. A non-Hermitian Hamilton operator and the physics of open quantum systems. J. 
Phys. A 42 , 153001 (2009). 

11. Heiss, W. The physics of exceptional points. J. Phys. A 45 , 444016 (2012). 

12. Bender, C. M. & Boettcher, S. Real spectra in non-Hermitian Hamiltonians having PT sym¬ 
metry. Phys. Rev. Lett. 80 , 5243 (1998). 

13. Klaiman, S., Gunther, U. & Moiseyev, N. Visualization of branch points in PT-symmetric 
waveguides. Phys. Rev. Lett. 101 , 080402 (2008). 


11 



14. Riiter, C. E. el al. Observation of parity-time symmetry in optics. Nature Physics 6, 192-195 
( 2010 ). 

15. Chong, Y., Ge, L. & Stone, A. D. PT-symmetry breaking and laser-absorber modes in optical 
scattering systems. Phys. Rev. Lett. 106 , 093902 (2011). 

16. Guo, A. et al. Observation of PT-symmetry breaking in complex optical potentials. Phys. Rev. 
Lett. 103 , 093902 (2009). 

17. Makris, K., El-Ganainy, R., Christodoulides, D. & Musslimani, Z. H. Beam dynamics in PT 
symmetric optical lattices. Phys. Rev. Lett. 100 , 103904 (2008). 

18. Longhi, S. Bloch oscillations in complex crystals with PT symmetry. Phys. Rev. Lett. 103 , 
123601 (2009). 

19. Lin, Z. et al. Unidirectional invisibility induced by PT-symmetric periodic structures. Phys. 
Rev. Lett. 106,213901 (2011). 

20. Regensburger, A. et al. Parity-time synthetic photonic lattices. Nature 488 , 167-171 (2012). 

21. Feng, L. et al. Experimental demonstration of a unidirectional reflectionless parity-time meta¬ 
material at optical frequencies. Nature materials 12 , 108-113 (2013). 

22. Peng, B. et al. Parity-time-symmetric whispering-gallery microcavities. Nature Physics 10 , 
394-398 (2014). 

23. Chang, L. et al. Parity-time symmetry and variable optical isolation in active-passive-coupled 
microresonators. Nature Photonics 8 , 524-529 (2014). 


12 



24. Liertzer, M. el al. Pump-induced exceptional points in lasers. Phys. Rev. Lett. 108 , 173901 
( 2012 ). 

25. Brandstetter, M. et al. Reversing the pump dependence of a laser at an exceptional point. 
Nature communications 5 (2014). 

26. Peng, B. et al. Loss-induced suppression and revival of lasing. Science 346, 328-332 (2014). 

27. Hodaei, H., Miri, M.-A., Heinrich, M., Christodoulides, D. N. & Khajavikhan, M. Parity¬ 
time-symmetric microring lasers. Science 346, 975-978 (2014). 

28. Feng, L., Wong, Z. J., Ma, R.-M., Wang, Y. & Zhang, X. Single-mode laser by parity-time 
symmetry breaking. Science 346, 972-975 (2014). 

29. Cao, H. & Wiersig, J. Dielectric microcavities: Model systems for wave chaos and non- 
Hermitian physics. Rev. Mod. Phys. 87 , 61 (2015). 

30. Joannopoulos, J. D., Johnson, S. G., Winn, J. N. & Meade, R. D. Photonic Crystals: Molding 
the Flow of Light (Princeton University Press, 2008), 2 edn. 

31. Sakoda, K. Proof of the universality of mode symmetries in creating photonic dirac cones. 
Opt. Express 20 , 25181-25194 (2012). 

32. Chan, C., Hang, Z. H. & Huang, X. Dirac dispersion in two-dimensional photonic crystals. 
Advances in OptoElectronics 2012 (2012). 

33. Mei, J., Wu, Y., Chan, C. & Zhang, Z.-Q. First-principles study of Dirac and Dirac-like cones 
in phononic and photonic crystals. Phys. Rev. B 86, 035141 (2012). 


13 



34. Lee, J. et al. Observation and differentiation of unique high-0 optical resonances near zero 
wave vector in macroscopic photonic crystal slabs. Phys. Rev. Lett. 109, 067401 (2012). 

35. Dembowski, C. et al. Experimental observation of the topological structure of exceptional 
points. Phys. Rev. Lett. 86, 787-790 (2001). 

36. Miroshnichenko, A. E., Flach, S. & Kivshar, Y. S. Fano resonances in nanoscale structures. 
Rev. Mod. Phys. 82, 2257 (2010). 

37. Hsu, C. W., DeLacy, B. G., Johnson, S. G., Joannopoulos, J. D. & Soljacic, M. Theoretical 
criteria for scattering dark states in nano structured particles. Nano Lett. 14, 2783-2788 (2014). 

38. Fleischhauer, M., Imamoglu, A. & Marangos, J. P. Electromagnetically induced transparency: 
Optics in coherent media. Rev. Mod. Phys. 77, 633 (2005). 

39. Hsu, C. W. et al. Observation of trapped light within the radiation continuum. Nature 499, 
188-191 (2013). 

40. Suh, W., Wang, Z. & Fan, S. Temporal coupled-mode theory and the presence of non- 
orthogonal modes in lossless multimode cavities. IEEE J. Quantum Electron. 40, 1511-1518 
(2004). 

41. Wiersig, J. Enhancing the sensitivity of frequency and energy splitting detection by using 
exceptional points: application to microcavity sensors for single-particle detection. Phys. Rev. 
Lett. 112, 203901 (2014). 


14 



42. Yoo, G., Sim, H.-S. & Schomerus, H. Quantum noise and mode nonorthogonality in non- 
Hermitian PT-symmetric optical resonators. Pliys. Rev. A 84 , 063833 (2011). 

43. Chua, S.-L., Lu, L., Bravo-Abad, J., Joannopoulos, J. D. & Soljacic, M. Larger-area single¬ 
mode photonic crystal surface-emitting lasers enabled by an accidental dirac point. Opt. Lett. 
39, 2072-2075 (2014). 

44. Szameit, A., Rechtsman, M. C., Bahat-Treidel, O. & Segev, M. PT-symmetry in honeycomb 
photonic lattices. Phys. Rev. A 84 , 021806 (2011). 

45. Lu, L., Joannopoulos, J. D. & Soljacic, M. Topological photonics. Nat Photon 8, 821-829 
(2014). 

46. Hsu, C. W. et al. Bloch surface eigenstates within the radiation continuum. Light: Science & 
Applications 2, e84 (2013). 

47. Haus, H. A. Waves and Fields in Optoelectronics (Prentice-Hall, Englewood Cliffs, NJ, 1984). 

48. Fan, S. & Joannopoulos, J. D. Analysis of guided resonances in photonic crystal slabs. Phys. 
Rev. B 65, 235112 (2002). 

49. Fan, S., Suh, W. & Joannopoulos, J. D. Temporal coupled-mode theory for the Fano resonance 
in optical resonators. J. Opt. Soc. Am. A 20, 569-572 (2003). 

Supplementary Information is available in the online version of the paper. 

Acknowledgments The authors thank Dr. Tim Savas for fabrication of the samples. Also, the authors 

thank Fan Wang, Yi Yang, Nick Rivera, Scott Skirlo, Dr. Owen Miller, and Prof. Steven G. Johnson for 


15 



helpful discussions. This work was partly supported by the Army Research Office through the Institute 


for Soldier Nanotechnologies under contract no. W911NF-07-D0004 and no. W911NF-13-D-0001. B.Z., 
L.L., and M.S. were partly supported by S3TEC, an Energy Frontier Research Center funded by the US 
Department of Energy under grant no. de-sc0001299 L.L. was supported in part by the Materials Research 
Science and Engineering Center of the National Science Foundation (award no. DMR-1419807). I.K. was 
supported in part by Marie Curie grant no. 328853-MC-BSiCS. 

Author Contributions All authors discussed the results and made critical contributions to the work. 

Author Information Reprints and permissions information is available at www.nature.com/reprints. The 
authors declare no competing financial interests. Correspondence and requests for materials should be 
addressed to B.Z. (email: bozhen@mit.edu). 

Competing financial interests The authors declare no competing financial interests. 


16 



Hermitian 



0.55 


a- 0.54 


0.01 


0 

\k\gJ2n 


0.01 


O 

K 

CM 

3 


CD 

=3 

cr 

CD 


b) 


non-Hermitian 

Re(co) c) 


lm(co) 


M -- r -- X 


M -- r -- X 





Figure 1: Accidental degeneracy in Hermitian and non-Hermitian photonic crystals (PhC). 

a, Band structure of a 2D PhC consisting of a square lattice of circular air holes. Tuning the 
radius r leads to accidental degeneracy between a non-degenerate quadrapole band and two doubly 
degenerate dipole bands, resulting in two bands with linear Dirac dispersion (red and blue) and a 
flat band (gray). b,c, The real and imaginary parts of the eigenvalues of an open, and therefore non- 
Hermitian, system: a PhC slab with finite thickness h. By tuning the radius, accidental degeneracy 
in the real part can be achieved, but the Dirac dispersion is deformed due to the non-Hermiticity. 
The analytic model predicts that the real (imaginary) part of the eigenvalue stays as a constant 
within (outside) a ring in the wavevector space, indicating two flat bands in dispersion, with a ring 
of exceptional points (EPs) where both the real and the imaginary parts are degenerate. In the upper 
panels, solid lines are from the analytic model and symbols are from numerical simulations: red 
squares represent the band connecting to the quadrapole mode at the center; blue circles represent 
the band connecting to the dipole mode at the center; and gray crosses represent the third band that 
is decoupled from the previous two due to symmetry. The 3D plots in the lower panels are from 
simulations. 
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Figure 2: Experimental reflectivity spectrum and accidental Dirac dispersion, a, SEM images 
of the PhC samples: side view (upper panel) and top view (lower panel), b, Schematic drawing 
of the measurement setup. Light from a super-continuum source reflects off the PhC slab and 
is collected using a spectrometer. The incident angle is controlled using a precision rotationary 
stage. (BS: beam splitter; SP: spectrometer) c, Reflectivity spectrum of the sample measured 
experimentally (upper panel) and calculated numerically (lower panel) along the T to X and the F 
to M directions. The peak location of reflectivity reveals the Hermitian part of the system, which 
forms Dirac dispersion due to accidental degeneracy. White lines in the lower panel indicate real 
part of the eigenvalues, d, Three line cuts of reflectivity from simulation results. Also shown are 
the complex eigenvalues (hollow circles) calculated numerically. At large angles (0.8°), the two 
resonances are far apart, so the reflectivity peaks (red arrows) are close to the actual positions of 
the complex eigenvalues. However, at small angles (0.3°, 0.1°), the coupling between resonances 
cause the resonance peaks (red arrows) to have much greater separations in frequencies compared 
to the complex eigenvalues. The black arrows mark the dips in reflectivity that correspond to the 
coupled-resonator induced transparency (CRIT, see text for details). 
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Figure 3: Experimental demonstration of exceptional ring, a, Examples of reflection spectrum 
from the sample at three different angles (0.8° blue, 0.3° green and 0.1° magenta, solid lines) 
measured with s-polarized light along F to X direction (same setup as in numerical simulations 
shown in Fig. 2d), fitted with the TCMT expression (equation ( |S.15| )) (black dashed lines). At 
each angle, the position of the complex eigenvalue extracted experimentally are shown as hollow 
circles, b, Complex eigenvalues extracted experimentally (symbols), with comparison to numerical 
simulation results (dashed lines) for both the real part (left panel) and the imaginary part (right 
panel). Red squares and dashed lines are used for the band with zero radiation loss at the T point, 
blue circles and dashed lines for the band with finite radiation loss at the T point, and gray crosses 
and dashed lines for the third band decoupled from the previous two due to symmetry, c, Positions 
of the eigenvalues approach and become very close to each other (indicated by the two brown 
arrows), demonstrating near EP features in different directions in the momentum space and the 
existence of an exceptional ring. 


19 




























Supplementary information 

Section I. Effective Hamiltonian of accidental Dirac points in Hermitian systems 

To the leading order of approximation, the effective Hamiltonian for accidental Dirac cones in 
Hermitian systems (2D PhC) is written as a 3 x 3 matrix due to the involvement of three bands: 


w 0 

VgK 

Vgky 

Vgk X 

Wo 

0 

Vgky 

0 

W 0 


that can be transformed into: 


UHgU 1 


with the orthogonal transformation matrix 


UJ Q Vgk 0 

Vgk w 0 0 
0 0 cjq 


1 0 0 

U — 0 cos 9 sin 9 (S.3) 

^ 0 — sin 9 cos 6 ^ 

Here, cos 6 = k x /k, sin 9 = k y /k, k x . y are in-plane wavevectors. After transformation, the 3 x 3 
matrix becomes two isolated blocks: the upper 2x2 block gives the conical Dirac dispersion 
(w = w 0 ± Vgk), while the lower block is the intersecting flat band (w = w 0 ). 

Section II. Effective non-Hermitian Hamiltonian of the exceptional ring 

For a 3D PhC slab that has finite thickness, the two dipole modes become resonances with finite 
lifetime due to their coupling to radiation; therefore, their eigenvalues become complex (w 0 — v'ov-i)• 
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With C 4 rotational symmetry, these two dipole modes are identical to each other under an 90 c 


rotation and therefore share the same complex eigenvalue. Meanwhile, the quadrapole mode does 
not couple to radiation at the T point due to symmetry mismatch, and to leading order its eigenvalue 
remains at cu 0 . The effective non-Hermitian Hamiltonian of the 3D PhC slab becomes 


/ 


rr3D _ 

-^eff — 


Vgky 


UJfl Vg k x 

Vgk x u 0 - *7d o 

0 Uo - i7d 




Vgky 


(S.4) 


which transforms to 


/ 


uh™u t = 


uj 0 v g k 

Vgk UJ 0 - ?7d 


V 


0 


0 


0 
0 

ujo - ^ 7 d 


\ 


(S.5) 


with the same matrix U as in equation S.3 The upper 2x2 block is the H cS we refer to in 
equation (|T|) that givies rise to an exceptional ring, while the lower block is the intersecting flat 
band. 

Section III, Existence of exceptional points along every direction in momentum space 

In this section, we demonstrate that EPs exist in all directions in the k space, not only for a sim¬ 
plified Hamiltonian (equation jT[), but also for realistic structures. To prove their existence, we use 
the unique topological property of EPs: when the system evolves adiabatically in the parameter 
space around an EP, the eigenvalues will switch their positions at the end of the loop 10 35 . In our 
system, the parameter space in which we choose to evolve the eigenfunctions is three-dimensional, 
consisting of the two in-plane wavevectors (k x , k y ) and the radius of the air holes r, as shown in 
Fig. S2a. Here, r can also be other parameters, like the refractive index of the PhC slab (n), the 
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periodicity of the square lattice (a), or the thickness of the slab (. h ). For simplicity of this demon¬ 
stration, we choose r as the varying parameter while keeping all other parameters (n, a, and h ) 
fixed throughout. 

First, we compare the evolution of the eigenvalues when the system parameters follow (1) a 
loop that does not enclose an EP, and (2) one that encloses an EP. Following the loop A —>• B —>■ 
C —> D —> A in Fig. S2a,b that does not enclose an EP (point E P ), we see that the complex 
eigenvalues come back to themselves at the end of the loop (Fig. S2c where the red dot and the 
blue dot return to their initial positions at the end of the loop). However, following the loop 
A' —>■ B' —>• C' —» D' —>• A! in Fig. S2d, which encloses an EP (point E P ), we see that the 
complex eigenvalues switch their positions in the complex plane (Fig. S2f where the red dot and 
the blue dot switch their positions). This switching of the eigenvalues shows the existence of an 
EP along the Y to X direction, at some particular value of k x and some particular value of radius 
r. This shows the existence of the EP without having to locate the exact parameters of k x and r at 
which the EP occurs. 

Similarly, we can evolve the parameters along any direction 9 = tan (k y /k x ) in the k space 
and check if an EP exists along this direction or not. As two examples, we show the evolution of 
the complex eigenvalues when we evolve the parameters along the 9 = n/8 direction following 
the loop A" —» B" —» C" —» D" —* A" and along the 9 — 7t/ 4 direction following the loop 
A"' —» B'" —> C" —> D"' —> A'" in Fig. S2g,h. In both cases, we observe the switching of the 
eigenvalues, showing the existence of an EP along these two directions. The same should hold for 
every direction in k space. 
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The above calculations show that for every direction 6 we examined in the k space, there 


is always a particular combination of k c and r c , which supports an EP. However, we note that in 
general, different directions can have different k c and different r c , so the exceptional ring for the 
realistic PhC slab structure is parameterized by k c {9 ) and r,,{()). This angular variation of k c (Q) 
and r c (9 ) can be described by introducing higher order corrections in the effective Hamiltonian, 
which we examine in the next section. 

Section IV, Generalization of the effective Hamiltonian 

Here, we generalize the effective Hamiltonian in equation (|TJ) and ( |S.5[ ). First, the radiation of the 
quadrapole mode is zero only at the T point; away from the T point, the quadrapole mode has a 
/. -dependent radiation that is small but non-zero, which we denote with y q . Second, we consider 
possible deviation from accidental degeneracy in the Hermitian part, with a frequency walk-off 5. 
With these two additional ingredients, the effective Hamiltonian becomes 


/ , 7 \ 

uj 0 + o Vgk 


V 


Vgk w 0 


( \ 

7q \/7q7d 

v x/T^Td 7d j 


(S.6) 


with complex eigenvalues of 

. 5 .7q + 7d 

= Wo + - - l —-— 


± \ {vgk-i^/% 7d) 2 - 


7d - 7q 


— I — 
2 


(S.7) 


which generalizes equations (|T) and (|2]). We note that the off-diagonal term x/7q7d in equa¬ 
tion (|S.6[) is required by energy conservation and time-reversal symmctr) i 4(J l 46 ! as we will discuss 


more in the next section. Equation ( |S.7[ ) shows that EP occurs when the two conditions 

k = (7d - 7q)/( 2 v g) ~ 7d/2 v g , 

& = 2 x/7dOh, 


(S.8) 
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are satisfied. In the region of momentum space of interest, 7 q is much smaller than 70 (this can 
be seen, for example, from the imaginary parts of Fig. Sla,c), so the first condition becomes 
hjQ « 7 d / 2 v g , same as in the simplified model. For a given direction 9 in the k space (as discussed 
in the previous section), we can vary the magnitude k and the radius r to find the k c (9 ) and r c {9 ) 
where these two conditions are met simultaneously. 

We can now analyze the angular dependence of k c {9 ) and r,,{()) without having to find their 
exact values. The first condition of equation ( |S. 8 [ ) says that the angular dependence of k c {9 ) comes 
from 7 d and v g ; in the PhC slab structure here, we find that 7 ,] varies by about 20% as the angle 
6 = tan (k y /k x ) is varied; while v g remains almost the same; therefore, k c {9 ) potentially varies by 
around 10% along the exceptional ring. For the second condition of equation ( |S. 8 [ ), we have 7,1 ~ 
5 x 10“ 3 o;o and 7 q « 5 x 10” 5 o;o for our PhC slab structure, so S c = 2 ^/ 777 ^ « lx 10” 3 o;o- Again, 
7 d and 7 q vary by about 20% as the angle 0 is changed, so S c can change by around 2 x 10 4 u 0 . 
Empirically, we find that a change of 5 by 2 x 10 ~ 4 cu'o corresponds to a change in the radius r of 
around 0.06 nm, which is the estimated range of variation for r c {9 ) of all 9 6 [0, 2n). This angular 
variation is much smaller than our structure can resolve in practice, since the radii of different 
holes within one fabricated PhC slab will already differ by more than 0.06 nm. So, in practice a 
given fabricated structure can be close to EP along all different directions 8, but is unlikely to be 
an exact EP for any direction. 

Section V, Temporal Coupled Mode Theory (TCMT) 

To connect the Hamiltonian of the resonances to the experimentally measured reflectivity, we resort 
to temporal coupled-mode theory (TCMT ) 30 47 . Here, we consider a very general setup with an 
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arbitrary number of resonances in the PhC slab. The time evolution of these n resonances, whose 
complex amplitudes are denoted by an n x 1 column vector A, is described by the Hamiltonian II 
and a driving term, 

A A 

— = -iHA + K V, (S.9) 

at 

where the Hamiltonian is an n x n non-Hermitian matrix 


H = Vt — iT — i7 nr , 


(S.10) 


with Q denoting its Hermitian part, — /T denoting its anti-Hermitian part from radiation loss, and 
—i 7 nr its anti-Hermitian part from non-radiative decays including absorption and surface rough¬ 
ness. For simplicity, we consider the same non-radiative loss for all resonances, so 7 nr is a real 
number instead of a matrix. 

Reflectivity measurements couple the n resonances to the incoming and outgoing planewaves, 
whose complex amplitudes we denote by two 2x1 column vectors, s + and s_. The direct reflection 
and transmission of the planewaves through the slab (in the absence of resonances) are described 
by a 2 x 2 complex symmetric matrix C, and 


s_ = Cs + + DA, 


(S.ll) 


where D and K in equation ( |S.9[ ) are 2 x n complex matrices denoting coupling between the 
resonances and the planewaves. We approximate the direct scattering matrix C by that of a ho¬ 
mogeneous slab whose permittivity is equal to the spatial average of the PhC slab ) 39 l 48 E^. Lastly, 
outgoing planewaves into the silica substrate are reflected at the silica-silicon interface, so 


s 2 + = e 2l ^ hs r 23 S 2 -, 


(S.12) 
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where h s is the thickness of the silica substrate with refractive index n s = 1.46, /3 = ri^uS 1 jc 1 — | k 
is the propagation constant in silica, and r 23 is the Fresnel reflection coefficient between silica and 
the underlying silicon. The formalism described above is the same as Ref. [39] except that here we 
describe the n resonances in a more general setting that accounts for their coupling (off-diagonal 
terms of H ) and therefore their non-orthogonality. 


For steady state with e ,aJ time dependence, we solve for vector A from equation (S.9) to 


get the scattering matrix of the whole system that includes both direct and resonant processes, 


S~ — (C + C res )>S+, 


(S-13) 


where the effect of the n resonances is captured in a 2 x 2 matrix 


C re s = lD(uj- H)~^ 


(S.14) 


We can solve equation (JSTT2J) and equation ( |S. 1 3[ ) to obtain the reflectivity 

2 


R 


TCMT 


Sl- 


« 1 + 


(S.15) 


In this expression, the only unknown is G' res . Therefore, by comparing the experimentally mea¬ 


sured reflectivity spectrum R(u) and the one given by TCMT in equation flS.15[ ), we can extract 
the unknown parameters in the resonant scattering matrix C res and obtain the eigenvalues of the 
Hamiltonian H. 

The remaining task is to write C res using as few unknowns as possible so that the eigenvalues 


of H can be extracted unambiguously. In equation (S.14), there are a large number of unknowns 


in the matrix elements of H, D, and K, but there is much redundancy because the matrix elements 
are not independent variables and because Cres is independent of the basis choice. Below, we show 
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that we can express C ies with only 2n + 1 unknown real numbers, and these 2 n + 1 real numbers 


are enough to determine the n complex eigenvalues of H. 

First, we normalize the amplitudes of A and s± such that their magnitude squared are the 
energy of the resonances per unit cell and the power of the incoming/outgoing planewaves per unit 
cell, respectively. Then, energy conservation, time-reversal symmetry, and C 2 rotational symmetry 
of the PhC slat 4 '® require the direct scattering matrix to satisfy = C* = C~ l and the coupling 
matrices to satisfy D' 1) = 2T, K = D, and CD* = —D. It follows that the matrix T is real 
and symmetric. Next, using the Woodbury matrix identity and these constrains, we can rewrite 


equation (S.14) as 


C res = -2 W (2 + W)~ 1 C, 


(S.16) 


where W = iD{u — Q + i^)~ l D^ is a 2-by-2 matrix. We note that the matrix W, and therefore 
the matrix C ies , is invariant under a change of basis for the resonances through any orthogonal 
matrix U (where Q is transformed to UilU~ x , and I) is transformed to DU~ l ). Therefore, we are 
free to choose any basis. Given the expression for W, we choose the basis where Q is diagonal, so 
Vtij = fljSij, with being the eigenvalues of fl. 

To proceed further, we note that the PhC slab sits on a silica substrate with 7i s = 1.46 and is 
immersed in a liquid with n = 1.48, so the structure is nearly symmetric in z direction. The mirror 
symmetry requires the coupling to the two sides to be symmetric or anti-symmetric 30 , 

^-=a j = ± 1, j = l (S.17) 

where a :l = 1 for TE-like resonances and a 3 = —1 for TM-like resonances, in the convention 
where (E x , E y ) determines the phase of A, and s±. Then, the diagonal elements of T are related 
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to D by Tjj = 7 j = \Dij\~, and in this basis we have 




* 7 j 


( \ 

1 o-j 


(S.18) 


\ a> V 


^ UJ - fL + i7n 

7=1 J 

This completes our derivation. Equations ( |S.16[ ) and ( |S. 1 8[ ) provide an expression for C res that 
depends only on 2 n + 1 unknown non-negative real numbers: the n eigenvalues {Q ? }” =] of the 
Hermitian matrix fl, the n diagonal elements {77 ) 7=1 of the real-symmetric radiation matrix T in 
the basis where is diagonal, and the non-radiative decay rate 7 nr . 

At each angle and each polarization, we fit the experimentally measured reflectivity spectrum 


R(u) to the TCMT expression equation ( |S. 1 5[ ) to determine these 2n + 1 unknown parameters. 
Fig. 3a and Fig. S 6 a show the comparison between the experimental reflectivity spectrum and 
the fitted TCMT reflectivity spectrum at some representative angles. The near-perfect agreement 
between the two demonstrates the validity of the TCMT model. 

To obtain the eigenvalues of the Hamiltonian H, we also need to know the off-diagonal 
elements of T. From I) ' I) = 2T and D \ 3 / D 2j = Oj , we see that Y, j — 0 when resonance % and 
resonance j have different symmetries in z (i.e. when 707 7 ^ 1 ), and that T I? = 7^/777 when 
77 = 1. In the latter case, the sign of T^ depends on the choice of basis; the eigenvalues of H 
are independent of the basis choice, so to calculate the eigenvalues of H, we can simply take the 
positive root for all of the non-zero off-diagonal elements of T. 

We note that the model Hamiltonians introduced previously, such as equation (|T) in the main 


text and equations (S.5) and (S. 6 ) above, are all special cases of the general Hamiltonian in equa¬ 


tion (S.10) that we consider in the TCMT formalism in this section. Those model Hamiltonians fix 


the number of resonances, assume simple forms of their parameters, and choose a specific basis 
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in order to convey the physical picture. Meanwhile, the TCMT formalism in this section does not 


make such assumptions (aside from basic principles such as energy conservation and time-reversal 
symmetry) so that it can be used as an unbiased method for analyzing the experimental data. 


We also note that the TCMT equations, from (S.9) to (S.17), are all written in the general 


matrix notation where one is free to choose any basis for the Hamiltonian II: we only make the 


specific basis choice (the basis where il is diagonal) in equation ( |S.18[ ) in order to simplify the 
expression for matrix W, and in equation ([3]) of the main text in order to emphasize the eigenvalues 
il. Meanwhile, physical observables, such as C res in equation < |S.14[ ), Rtcmt in equation ( |S.15[ ), 
and the eigenvalues of H, are all independent of the basis choice. 

Section VI, Reflection peaks and CRIT 

In this section, we use a simplified scenario (a special case of the previous section) to illustrate 
that the peaks of the reflectivity generally follow the eigenvalues of f) and to show the coupled- 
resonator-induced-transparency (CRIT). 

Consider a simplified scenario with two resonances of the same symmetry in z and without 
non-radiative loss (i.e. n — 2, a\ — a 2 , 7 „ v = 0), and ignore the direct Fresnel reflection between 
the dielectric layers (so that the direct scattering matrix C has no reflection, and that s 2 + = 0 in 
equation ( |S.12[ )). In such case, equations ( |S.15[ ) ( |S.16[ ) ( |S.18[ ) give 


-Rtcmt(^) - 1 f0( v “77 \ — —V ^ -(S- 19 ) 

1 + /-(tu) j[CO) LO — ill UJ — il 2 

We immediately see that the reflectivity reaches its maximal value of 1 when co = Hi or u = il 2 , 
namely at the eigenvalues of the matrix il. Another feature we can observe is that the reflectivity is 
0 when u> = ( 71 H 2 + 72 ^ 1)/(71 + 72 ), which is a phenomenon called coupled-resonator-induced- 
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transparency (CRITpM 37 !. 

We emphasize that the reflectivity peaks are different from the complex eigenvalues of the 
Hamiltonian H. Consider a simple example with fi 1;2 = cu 0 ± b, and 71 = 72 = b. The reflection 
peaks at f ) lj2 = u 0 ±b, while the two complex eigenvalues are degenerate at uj + = cu_ = cu 0 — ib, 
whose real part is in the middle of the two reflection peaks. This explains the reflectivity from the 
PhC slabs at 0.3° shown in Fig. 2d and Fig. 3a, where the degenerate complex eigenvalues of the 
system are in between the two reflection peaks. 

In Fig. S5, we use some examples to illustrate the difference between the reflectivity peaks 
and the complex eigenvalues. Fig. S5a shows the case when there is only one resonance (removing 


one of the two terms in equation (S.19)) with complex eigenvalue uq— * 7 ; in this case, the reflection 


peak (red arrow) is at the same position as the real part of the complex eigenvalue. In contrast, 
Fig. S5b shows the case when there are two resonances (equation ( |S.19[ )) with f)i j2 fixed at cu 0 ± b; 
as we vary 71 , 2 , the complex eigenvalues (circles) vary accordingly, whereas the reflectivity peaks 
(red arrows) always show up at fi 12 . 

For the realistic PhC slab structure in our experiment, the reflectivity is described by the more 


general expression, equations (S.15), but equation (S.19) serves as a qualitative approximation near 


the frequency range of interest, because the far-away resonances do not contribute much, the non- 
radiative loss is small, and the Fresnel reflection between the dielectric layers (liquid, Si 3 N 4 , silica, 
and silicon) is small. So, we can still see the general trend that the reflectivity peaks follow the 
eigenvalues of Q (as evident by comparing Fig. S3 and Fig. S 6 b), and we can still see reflectivity 
dips for CRIT (such as in Fig. 3a). 


30 










a) r < r c 

m -— r —- k 




FIG. SI: Simulation results of the complex eigenvalues of the PhC slabs with and without 
accidental degeneracy. The real (upper panels) and imaginary (lower panels) parts of the complex 
eigenvalues are shown for structures with (b) and without (a, c) accidental degeneracy. The bands 
with quadrapole modes in the middle of the Brillouin zone are shown in red solid lines, while the 
bands with the dipole mode are shown in blue solid lines and gray dashed lines. The bands shown 
in red and blue solid lines couple to each other, while the band in gray dashed lines is decoupled 
from the other two due to symmetry. When accidental degeneracy happens (at r = r c as shown in 
b), the characteristic branching features are observed demonstrating the existence of EPs. When 
the accidental degeneracy is lifted, the quadrupole band splits from the dipole bands from different 
directions: from the bottom when radius of the holes is too small (r < r c as shown in a), or from 
the top when the radius is too big (r > r c as shown in c). 
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Eigenvalues stay when not looping around EPs 
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Eigenvalues switch when looping around EPs 
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FIG. S2. Existence of EP along every direction in the momentum space for the realistic 
PhC slab structure, a, A loop is created in the parameter space of the structure (A —» B —>• C —>• 
D —>• A), which does not enclose the EP of the system (point E p ). Here, r is the radius of the 
air holes, and k XjV are the in-plane wavevectors. b, For each point along the loop, we numerically 
calculate the eigenvalues of the PhC slab with the corresponding hole radius at the corresponding 
in-plane wavevector. c, The complex eigenvalues return to their initial positions at the end of the 
loop (namely, the blue dot and the red dot come back to themselves) when the system parameters 
come back to point A. d,e,f, Another loop is created ( A' —tB'—tC'—tD'—} A'), which 
encloses an EP of the system (the same point E p as in a.). Following this new loop, the two 
eigenvalues switch their positions at the end of the loop (namely, the blue dot and the red dot 
switch their positions) when the system parameters come back to point A! . g,h. The two complex 
eigenvalues always switch their positions when we choose the right loops along other directions in 
the momentum space (T to N in g and F to M in hj. 


33 



Reflectivity spectrum 



FIG. S3: Experimental results of reflectivity showing an accidental Dirac cone. Light with 
different polarizations (s and p) is selected to excite different resonances of the PhC slab along 
different directions in the k space. Depending on the choice of polarization, the two bands forming 
the conical dispersion are excited (s-polarized along F-X and p-polarized along F-M, shown in the 
left panel), or the flat band in the middle is excited (p-polarized along F-X and s-polarized along 
T-M, shown in the right panel). 
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FIG. S4: Experimental results of reflectivity from PhC slabs with and without accidental 
degeneracy. Angle-resolved reflectivity along the T to X direction, measured for three different 
PhC slabs: a, with smaller hole radius than the structure with accidental degeneracy (r < r')\ b, 
with accidental degeneracy (r = r'); c, with bigger hole radius than the structure with accidental 
degeneracy (r > r'). The reflectivity peaks of the structures without accidental degeneracy (a,c) 
follow quadratic dispersions; while the reflectivity peaks of the structure with accidental degen¬ 
eracy (b) follow linear Dirac dispersion. Data shown in (b) is the same data as in Fig. 2c the left 
panel of Fig. S3. 
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a) Reflection of one resonance 



b) Reflection of two coupled resonances 



FIG. S5: Illustrative reflectivity spectrum from one resonance and from two coupled reso¬ 
nances. a, When a single resonance dominates, the reflectivity peak is at the same position as the 
real part of the complex eigenvalue, b, With two coupled resonances, the reflectivity peaks (red 
arrows) no longer follow the eigenvalues of the system (blue circles). As we vary the radiation loss 
7 i ) 2 of the two resonances while fixing the eigenvalues 0 12 of the Hermitian matrix, we see the 
complex eigenvalues vary whereas the reflectivity peaks are fixed. The middle panel of (b) shows 
a situation when the two complex eigenvalues coalesce into an EP. 
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b) Fitting results before diagonalization 
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FIG. S 6 : TCMT fitting and visualization of accidental Dirac cone, a, Examples of reflection 
spectrum measured at five different incident angles (0°, 0.14°, 0.24°, 0.5° and 1°) along the T- 


X direction for s polarization, with comparison to the TCMT expression in equation ( |S.15[ ) after 
fitting. Dotted lines indicate the resonances, with the two relevant resonances marked in red and 
blue, b, Parameters obtained from the TCMT fitting. The eigenvalues for the Hermitian part of 
the Hamiltonian, 0 1)2 , are shown in the left panel and reveal the Dirac dispersion arising from 
accidental degeneracy. The diagonal terms for the anti-Hermitian part of the Hamiltonian, 712 , 
are shown in the right panel. Note that the anti-Hermitian part of the Hamiltonian also has off- 
diagonal terms, so Oi )2 and 71 , 2 are not the eigenvalues of the Hamiltonian. The eigenvalues of the 
Hamiltonian are shown in Fig. 3 of the main text. 
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